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Abstract. By using a simple observation that the density pro- 
cesses appearing in Ito's martingale representation theorem are 
invariant under the change of measures, we establish a non-linear 
version of the Cameron-Martin formula for solutions of a class of 
systems of quasi-linear parabolic equations with non-linear terms 
of quadratic growth. We also construct a local stochastic flow 
and establish a Bismut type formula for such system of quasi- 
linear PDEs. Gradient estimates are obtained in terms of the 
probability representation of the solution. Another interesting 
aspect indicated in the paper is the connection between the non- 
linear Cameron-Martin formula and a class of forward-backward 
stochastic differential equations (FBSDEs). 
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1 Introduction 

The goal of the article is to present a unified approach for Girsanov's tech- 
niques of changing probability measures used in the recent literature on back- 
ward stochastic differential equations (BSDE). The framework we formulate 
has an advantage allowing us to bring together seemingly different sorts of re- 
sults on BSDE, forward-backward stochastic differential equations (FBSDE), 
and function equations on a probability space. 

Consider the following system of quasi-linear parabolic equations 



i = 1, • • • ,m, where the drift vector field involves a solution u and its to- 
tal derivative Vw = (f^j) through / = (P)j<d- In the interesting cases, 
/ J : M. d x ]R mxd M. (j = 1, • • • ,d) are Lipschitz continuous but unbounded, 
and therefore the non-linear term (often called the convection term) appear- 
ing in ( II. lj) is of quadratic growth. This is a kind of non-linear feature which 
appears in many physical PDEs, see for example [5] and [TT]. If m = d = 3 
and / J (w) = vP , then (jl.ip is a modification of the Navier-Stokes equations 
with the pressure term and the divergence-free condition dropped altogether, 
while the same non-linear convection term has been retained. This kind of 
PDEs has been used as simplified models for phenomena such as turbulence 
flows. Due to the special structure of the system (II. ID . the maximum princi- 
ple applies to \u(x, t)\ 2 , thus a bounded solution exists as long as the initial 
data u(x, 0) is regular and bounded, which makes a distinctive difference 
from the Navier-Stokes equations. According to Theorem 7.1 on page 596 in 
[H], if the initial data Uq is smooth and bounded with bounded derivatives, 
then a bounded, smooth solution u to the initial value problem of (II. ip ex- 
ists for all time. Our main interest in this article is to establish probabilistic 
representations for the solution u by applying Girsanov's theorem to BSDEs. 

To this end, it will be a good idea to look at Peng's non-linear Feymann- 
Kac formula (see [15]) for general quasi-linear parabolic equations with Lip- 
schitz non-linear term. The main idea in [TJ] and [15] can be described as 
following. Thanks to Ito's calculus, Brownian motion B = (B 1 , • ■ ■ , _B d ) may 
be considered as "coordinates" on the space of continuous paths equipped 
with the Wiener measure, and it is a much tested idea that one may "read 
out" solutions of quasi-linear partial differential equations in terms of B. 
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Suppose u(x,t) = ( , u 1 (x,t), • • • ,u m (x,t)) is a sufficiently smooth solution to 
the Cauchy problem of the system of quasi-linear parabolic equations 

-A -^-]u i + h*(u, Vu) = in R d x [0, oo) 
2 at J 



with the initial data u l (x, 0) = u l (x). Applying Ito's formula to Y t l = u l (x + 

dxi 

d r x 



B t , T — t) (and set Z\' j = ^j(x + B t , T - t)) to obtain 



Y T -Y} = E / + B„T - s)dB> 

= T f Z^dBl- [ T h\Y s ,Z s )ds 
i=i Jt Jt 



for t e [0,T]. The pair (Y, Z) is a solution to the stochastic differential 
equations 

d 

dYl = -h\Y t , Z t )dt + E z t j dB{ (1.2) 

i=i 

on (Q, J 7 , Tt-, P) with the terminal data Y T = u l (x + Bt). 

Peng [15] and Pardoux-Peng [TJ] made two crucial observations: firstly 
the process Z t may be recovered from the Ito representation of the martingale 
SI = Sq + f Y^j=i Z^dBi. Secondly 5* is indeed the martingale part of the 
continuous semimartingale Y, therefore (11.2j) is a closed system and may be 
solved by specifying a terminal value Yp. It was proved in Pardoux-Peng [TJ] 
that if h % are global Lipschitz continuous, then (11. 2p may be solved as long 
asY T E L 2 (tt,F T ,P). 

In the case that h l (y, z) has a special form such as P(Vi z)z^ , one is 

able to solve by firstly solve the trivial BSDE dYj = Y? j= i z t jdB t then 
change the probability measure, even for non-Lipschitz non-linear term h as 
long as / is Lipschitz continuous. We are thus able to extend the Cameron- 
Martin formula to a class of systems of quasi-linear parabolic equations with 
quadratic growth. The non-linear version of the Cameron-Martin formula, 
which may be considered as our contribution of this article, is of independent 
interest. 
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Suppose that the initial data u is Lipschitz continuous and bounded. Let 
(J~t)t>o be Brownian filtration, i.e. the filtration over W d of all continuous 
paths in R d generated by the coordinate process {B t : t > 0}, augmented by 
the Wiener measure. Let T > and x G R d be fixed but arbitrary. For each 
f = (f !, • ■ ■ , £ m ) where C G L°°(W d , ? T , P x ). Define 

B(Ot = B t - f f(Y(0s,Z(0s)ds 
Jo 

where Y(£) t = -E ,p (^|J r i) are bounded martingales and Z(£) = are 
determined by Ito's martingale representation: 

e-E p (e\Fo) = it f T z(a' j dB{. 

3=1 J ° 

We prove that there is a unique f G Jt, P 35 ) such that ^ = Uo(-B(f)r), 

and u(x,T) = E P (Y(^) ), which can be considered as the non-linear version 
of the Cameron-Martin formula. More information about the solution u may 
be obtained as we can represent the derivative J^- in terms of the process 
[Z^Y^t^Ti see Theorem 13.11 below. The non-linear version of the Cameron- 
Martin formula may be reformulated in terms of FBSDE as well, see (12.41) 
and Corollary 12.41 

The main reason we are interested in representations of solutions to phys- 
ical PDEs in terms of Brownian motion or in general in terms of Ito's diffu- 
sions, lies in the fact that it is then possible to employ probabilistic methods 
such as Ito's calculus, Malliavin's calculus of variations and path integra- 
tion method to the study of non-linear PDEs. We demonstrate this point 
by deriving explicit gradient estimates for solutions of a class of systems of 
quasi-linear parabolic equations. 

Finally we would like to point out that the type of BSDEs such as fll.2p 
with quadratic growth non-linear terms has been well studied in Kobylanski 
[7]. Her work has been extended and generalized substantially in Briand 
and Hu [3] and [2J. BSDEs with quadratic growth driven by martingales 
are solved by Morlais [13] and Tevzadze [TB] . However their methods can be 
applied to scalar BSDEs, not systems in general. 
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2 Girsanov's theorem, martingale represen- 
tation and BSDE 



Let B = (B 1 , ■ ■ • , B d ) be a Brownian motion in R d started at on a com- 
plete probability space (f^J 7 , P), (J~t)t>o be the filtration generated by B 
augmented by probability zero sets in T and = o{T t : t > 0}. Since 
we are interested in J r 00 -measurable random variables only, for simplicity, we 
assume that T = J-^. According to Ito's martingale representation theorem, 
any martingale S on (fl, J 7 , Tt, P) is continuous and has a unique represen- 
tation in terms of Ito's integration 



where Db(SY are predictable processes, called the density processes of S with 
respect to Brownian motion B. Since we will deal with several equivalent 
measures on (Q, J 7 ) at the same time, it is desirable to have labels associated 
with notations which involve a probability measure. We will follow this 
convention if confusions may arise. Therefore, E p and E p {-(J 7 *} denote the 
expectation and conditional expectation with respect to P respectively. 

Let Q be a probability measure on (O, J 7 ) equivalent to P, whose Radon- 
Nikodym's derivative with respect to P restricted on T% is denoted by R t , 
that is, | — R t f or t > q Then R is a positive martingale on (Q, J 7 , T t , P) 

with E p (R t ) = 1. Girsanov's theorem [6] establishes a correspondence be- 
tween local P-martingales and local Q-martingales. If X = (X t ) t > is a local 
martingale under probability P, then X is a continuous semi-martingale un- 
der Q and X t = X t - f* -^d(X } R) s is a local martingale under Q, where the 
bracket process (X, R) t is defined under the probability P, which is however 
invariant under the change of equivalent probability measure. 

It is convenient to formulate the Girsanov's transform in terms of expo- 
nential martingales. The exponential martingale S(N) t = expjiVt — h(N)t} 
of a continuous local martingale N (under P) is the unique solution to the 
stochastic exponential equation 



■Jo 

Up to an initial data R , the Radon-Nikodym derivative R = (Rt)t>o of Q 
with respect to the measure P must be an exponential martingale of some 
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continuous local martingale N, so that R t = Roexp{N t — ~(N) t }. Then 
X = X — (X, N) for any P-local martingale X. 

According to Levy's theorem, B = (B 1 , ■ ■ ■ ,B d ) is a Brownian motion 
under Q which is (J-" t )t>o-adapted, but the natural filtration (Ft)t>o of B may 
not coincide with (J r t)t>o- It can happen that T% is strictly smaller than J-j for 
some t. But, nevertheless any (J-" t )-martingale under Q may be represented 
as Ito's integral of (J-j)-predictable processes against B. The starting point 
of our approach is the following elementary fact. 

Lemma 2.1 For any local P '-martingale X (with respect to the filtration 
(Ft)t>o) Db{X) = Dg(X). That is, the density process appearing in Ito's 
martingale representation is invariant under the change of equivalent proba- 
bility measures. 

The lemma follows directly from the definition of density processes and 
the Girsanov's theorem, so its proof is left to the reader. With the help of 
this simple fact, we can deal with the question of how a backward stochastic 
differential equation is transformed under change of equivalent probability 
measures. 

Theorem 2.2 Suppose that g is global Lipschitz continuous, £ G L 2 (Q, Tt, P), 
and (Y, S) is the unique solution pair to BSDE 

dY t = -g(t, Y t , D B (S) t )dt + dS t , Y T = i (2.1) 

on (fi, J 7 , J*t, P) . LetS = S- (S, N) . Then ( Y, S) solves 

d 

dY t = -git, Y t , D s (S) t )dt + Dz(S) J t D B (N)idt + dS t (2.2) 

3=1 

on $l,JF,F t ,Q), Y T = i. 

Proof. Since (Y, S) is a solution to ( 12. ip . Y % are continuous semimartin- 
gale with martingale parts S j . By Lemma I2~T| D B (S l ) = D^(S l ), so that 

§ i - So = E fDzi&yM = J2 f D B {syM 

3=1 J ° 3=1 Jo 
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and 



Therefore 



3=1 J ° 



Y^-Yl = -J g\s,Y s ,D B (S) s )ds 

+ E f D^Sy s D B {N){ds + S T - § (2.3) 
i=i Ji 

which is valid under the probability Q. That is, the pair (Y, 5) solves BSDE 
( P on (fi, J 7 , Q), with terminal values K£ = f . ■ 

The most interesting case is the following special choice of N. 

Corollary 2.3 Under the same assumptions as in the previous theorem, and 

N t = Yl f P(Y s ,D B (S) s )dBi 
3=1 Jo 

where f : R m x R md R d is Borel measurable. Then (Y, S) solves the BSDE 



d 

dY t = -g(t, Y t , D s {3) t )dt + £ D B (S)lf(Y t , D B (S) t )dt + dS t 

3=1 

with the terminal value Yt = £ on (fl, J 7 , J-'t, Q)- 

One may reformulate Corollary 12. 3l in terms of forward-backward stochas- 
tic differential equations. Observe that, with the choice of N made in Corol- 
lary 12.31 X = x + B is the solution to the stochastic differential equation 

dX t = -f(Y s , Z t )dt + dB u X t = x 

while (Y, Z = D B (S)) is the solution of the BSDE 

dY t = -g(t, Y u Z t )dt + Z t dB u Y T = u (X T ), 
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thus (X, Y, Z) is a solution to the following forward-backward stochastic dif- 
ferential equations (FBSDEs) 



dX t = -f(Y s , Z t )dt + dB t , 
dY t = -g{t,Y t ,Z t )dt + Z t dB t , 
X t = x,Y T = u (X T ). 



(2.4) 



FBSDEs such as (12.41) have been studied by various authors, and are 
well presented in the research monograph [10J. By utilizing the fundamental 
apriori estimates established in [8], it has been proved in [5] that FBSDE 
(12.41) has a unique solution such that Y is bounded, as long as u is bounded 
and Lipschitz, and if / and g are Lipschitz continuous. 

Corollary 2.4 Let f and g be Lipschitz continuous, uq be bounded and Lip- 
schitz continuous. Let (X, Y, Z) be the unique solution of \2.J$ such that Y 
is bounded. Define Q on (Q, J 7 ?) by 



dQ 
IP 



expj^ f{Y s ,Z s ).dB s - l -J^ \f(Y s ,Z s )\ 2 d s y 



Then (Y,Z) is the unique solution (such that Y is bounded) to 



-g{t, Y u Z t ) + f 1 ^ ~ Z t)~ Z l ^ + Eli Z{dB{, 




on (O, J 7 , 7 U Q) 



(2.5) 



3 Cameron-Martin formula for PDEs 

In this section we apply the machinery developed in the previous section to 
the initial value problem of the following quasi-linear parabolic system 

— + Y,r(u,Vu)— = -Au> + g\u,Vu) inR d (3.1) 
i=i 

for i = 1, ■ ■ • , m, where / J and g l are global Lipschitz continuous in their 
arguments, together with the initial data u(x, •) = u (x) which is bounded 
and global Lipschitz continuous. According to Theorem 7.1 in [8J, there is 
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a unique bounded solution u(x, t) to the initial value problem of (13. ip . By 
Ito's formula, Y t = u(x + B t , T — t) and Z t = Vm(i + B t , T — t) solve the 
BSDE 



dYi 



-g\Y t ,Z t ) + J2f 3 (YuZt)Z^ 

3=1 



dt ■ £ Z'-dBi (3.2) 

3=1 



with the terminal value Y? = Uq(x + Bt) over the filtered probability space 
(fi, J 7 , J-i, P), which is indeed the unique bounded solution of BSDE (13.21) . 

On the other hand, since g l are global Lipschitz, according to Pardoux 
and Peng [14] . for any £ G L 2 (fi, J 7 ^, P) there is a unique solution pair (Y, Z) 
to the backward stochastic differential equation: 



dY t = -g(Y t , Z t )dt + zJ t dB l y t = Z 



(3.3) 



3=1 



on (Q, J 7 , Tti P), which is the system derived from (I3.2p with the quadratic 
non-linear term dropped . The solution (Y, Z) depends on the terminal value 
£, so it is denoted by Y(£) and Z(£) respectively. For x £ R d and T > 0, let 



and 



N(0t = it f 'f(Y(0s,Z(0s)dBi 
3=1 Jo 

B(0t = B t - I f(Y(0s,Z(0 s )ds. 
Jo 



(3.4) 



(3.5) 



Theorem 3.1 Let u be a classical solution of h3. 1\) on [0, T] x R d such that 
u(0, •) = uq(-). Suppose that £ e L°°(Q, Tt, P) is the solution to the function 
equation 

£ = u (x + B(Ot), 

then 

u(x + B(0t,T-t) = Y(0t 

and 

Vu(x + B(0t,T-t) = Z(0t 
for all t < T almost surely. 
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Proof. Define an equivalent probability measure Q by jp\ Tt = £(iV(£)) t . 
Then B(£) is a Brownian motion (up to time T) under Q. Let S = S(£) — 
(N(g),S(g)). According to Corollary E31 (Y(£),S) solves the following 
BSDE 

d 

dY t = Db{S)IP{Y u D s (S) t )dt - g(t, Y t , D § (S) t )dt + dS t 

3=1 

with terminal value Yt = <p(x + B(£)t) on (Q, J 7 , J-%, Q). This system is 
exactly the BSDE that Y t = u(x + B(£) t ,T — t) should satisfy, and the 
conclusions follow immediately. ■ 

By Theorem I3.1[ in order to provide a probabilistic representation for 
(I3.ip . the problem is reduced to solve the function equation 

Z = u (x + B T - J Z(£) s )ds\ := 0(0- (3.6) 

The following local existence of the solutions to (13. 6 j) is elementary. 

Proposition 3.2 Let C UQ and Cf be the Lipschitz constants for uq and f 
respectively, and 



T 



&c 2 uo c 2 



1 - 1. 



IfT<r, then 113.6}) admits a unique fixed point £ G L°°(Q, J-^, P). 
Proof. For £ and r] in L 2 (il, T T , P), 

10(0-0(^)1 <a„ f \f(Y(0s,Z(0s) ~ f(Y( V ) s ,Z( V ) s )\ds 
Jo 

<C U0 C f f \Y{C-v)s\ + m- V )s\ds 
Jo 

and furthermore, 

E[m)-m\ 2 ] <ci cjE 
< 



J \Y(Z- V ) s \ + \Z(Z-r)) s \ds\ 

2C 2 Uo C)T \ f{E\Y{i - V ) s \ 2 + E\Z{a - (3) s \ 2 )ds 
Uo 
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Now, for any 77 G L 2 (f2, jF t , P) 



\Y( v ) s \\l = E[\E(ri\F s )\ 2 } < IM 



2" 



and 



|Z(r/) s | 2 ds 



= ^[|y(^-y(77)or] 
<2|HI1, 



so that 



# [10(0 - m\ 2 ] < 2c 2 U0 c 2 f T \ f T ne - 1711^ + 2|ie - 



2^C 2 T(T + 2)||e-^|| 2 . 



Hence 



1 10(0 - <KV) 1 1 2 < V2C UQ C f VT(T + 2)||e-77|| 2 
and the claim follows from a simple application of the fixed point theorem. 



Theorem 3.3 Under the previous assumptions, the function equation ( Iff. 6]) 
has a unique solution. 

Proof. We try to extend the local solution constructed in the previous 
proposition to the case T > r. To do this we need a uniform gradient estimate 
of the solutions to PDE (13. ip : there exists a constant C u such that 

sup | Viz | (a;, t) < C u . 

(t,x)G[Q,T]xR d 

For the proof see for example [8] and [3]. Based on such constant C u we 
define 



1 



8C*CJ 



1 - 1. 



Now we construct the random variable £' on the time interval [r,r + r']. 
Consider the following function equation: 

£ = u(x + B' T+T , - £ +T f(s, Y(0s, Z(0s)ds, t j := 0'(O- (3.7) 
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where B' is a Brownian motion on [r, r + r'j denned by B s = B s — B T for 
s G [r,r + t']. Analogous to the proof of Proposition 13.21 (13. Tj) admits a 
unique fixed point £' G L°°(0, J> +T ', P) C L 2 (f2, Jv +r /, P). Based on such 
we have the following representation formulae on [r, r + t']: 

M (x + J B(Ot,2r + r'-t) = r(Ot 

and 

V«(a; + S(Ot,2r + T / -t) = Z(Ot 
for i G [t, r + t'], where 

s(0< = B 't- J* ^(0«)d«. 

We then move to the next interval [r + r', r + 2r'] and repeat the above 
procedure until we touch T. ■ 

4 Some applications 

In this section we establish some explicit gradient estimates for the solution of 
(13. ip by using the representation theorem 13.11 to demonstrate the usefulness 
of non-linear Cameron-Martin's formula. Further applications will appear in 
a separate paper. 

Let us retain the notations and assumptions established in the previous 
section. Let T > be fixed, u = (u 1 , • ■ ■ , u m ) be the unique solution to the 
initial value problem (13 .ip with initial data uq, where / (which determines 
the nature of the quadratic non-linear term) depends only on y, i.e. f(y, z) 
does not depend on z and g = 0. That is, u is the solution to the initial value 
problem of the following system of quasi-linear parabolic equations 

^ + y/i( u )*C = I A «* inR d (4.1) 

i=i 

for i — 1, • • • , m. 

Let £ be the solution to the function equation (13.61) established in Theorem 

Let Q be the equivalent measure with density process ■gp\ :F — R t — 
£(N(0h where 

N(0t = [ f(Y(0) s ).dB s . 



12 



Theorem 4.1 1) Let p G [1,2) and (Pt)t>o the heat semi-group, i.e. P t 
e^ tA . Then 



1^ P s \Vu i \ p (x,T-s)ds<d 1 -2 (E p (Y{O i )t)" Qf E p Rt*ds^J 



1-2 



for any t < T, i = 1, • • • , m. 
2) We have 



\Vuf (x,T) <\im tl0 ^E p (Y^) l )t 



for any t < T, % = 1, • • • , m. 
Proof. By Theorem 13.11 



so that 



(^(0*)*= / |V«f (x + B(0»,T-s)ds 
Jo 



and therefore 



E H I \Vu t \*(x + B(Os,t-s)ds = E 1J (Y(O i )f 
Jo 

On the other hand, for 1 < p < 2 we have 

EQ {Io \ Vui \ P ( x + B(£)s,T-s)ds^ 
= J E^{\Vuf(x + B(Os,T-s)}ds 

< jf (e q |i- 1 v«f (x + £(£)„, r - 2 (e q (i^)) 1 " 5 <*» 
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Since B(£) is a Brownian motion under Q, so that 

E Q [ \Vu i \ P {x + B{Os,T-s)ds 
Jo 

= [ P s (\Vuf(;T- S ))(x)ds 

Jo 

£ p 

< (E P \W\ 2 (:r + 5(O s ,T-s)) § ^E p Rf^j 2 ds 

< d 1 '2^E p \Vuf(x + B(0s,T-s)dsy^E p Rf~ p ds^j 

= d^{E p (Y{t)%f^E p Rf*d^j 2 

which is the first estimate. To prove the second one, we write the previous 
estimate as 

ft 



1 r 

- / P S |W| P (x,T - s)ds 
t Jo 



t 2 \l-f 



o 



Letting t — >■ one obtains that 

^J\vuf(x,T) < d^m t ^\J\E p (Y{tiy)t. 

then letting p f 2 we obtain 2). ■ 
Lemma 4.2 TTien /or any p G [1, 2) 



^(^')<exp^^ itaS J/fe)P 



Proof. Let 

H t = exp 



(2 
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which is exponential martingale, so that E p (H t ) = 1. Then 

r-t 



Rl~ v = H t exp 
< H t exp 

which yields the claim. ■ 
Corollary 4.3 We have 

P s IVtff (x,T - s)ds < exp 



P 



\f(Y(Os)\ 2 ds 



(2 - P? Jo 

— -t max \f(y)\ 2 

(2-p) 2 |y[<|«o|oo 



P 



2(2 -p) \y\<\u \ 



t max \f(y)\' 



(4.2) 

for any p G [1, 2), and i = 1, • • ■ , m, where P s = ea sA the heat semigroup. 
Proof. Observe that 

E P (Y{Q)s = E p [E p {Cm-E p {Q\F Q )\ 2 

= E p [E p (^ s y-E p (^ y] 

< I Kill 2 

^ l w olc» 

and therefore, the item 1) in Theorem 14.11 together with Lemma [4.21 yields 
the gradient estimate (14.21) . ■ 



5 Non-linear flow associated with quasi-linear 
PDEs 

In this section we construct a non-linear stochastic flow associated with the 
quasi-linear system (14 .11) . 

Assume that f j : R m ->• R d (j = 1, ■ ■ • , d) are differentiable with bounded 
1st, 2nd and 3rd derivatives, and the initial value u l (x,0) = u (x) {i = 
1, ■ • • ,m) are bounded, differentiable with 1st and 2nd bounded derivatives: 

|V fc-1 u | < C , |V fc /| < Ci for k = 1,2,3 

for some non-negative constants Co and C\, here V fc denote the fc-th deriva- 
tive in space variables. Of course V°wo = uo- I n general, we apply V to mean 
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the total derivative operator in space variables. For example, Vw means (f^-) 
but does not include the derivative in time parameter t. 

We are going to construct a continuous adapted process £ = (£ f ) with 
values in the function space C 2 (R d ,R m ) such that V j u(-,t) = E p (V j £ t ) at 
least for small t. The spirit in devising such a formula is quite similar to 
those initiated in the seminal works Bismut [T] and Malliavin [T2] . 

Consider the function space C([0, T]; C 2 (R d , R m ))). If 

£ G L°°(fi;C([0,!Zl;C£(R d ,R m ))) 

then for any w, G C([0, T]; C 2 (R d , R m )), so that £ t (u) G C 2 (R d ,R m ) 

and t -)■ £ t (u;) is continuous from [0,T] to C 2 (R d , R m ), and x -> 
has continuous and bounded first and second derivatives. Let H T denote 
the space of all functions £ in C([0, T]; C 2 (R d , R m ))) such that £ = 

(Ct)te[o,r] is adapted, iir is equipped with the L°°-norm, namely 

2 

||f|| =esssup sup SU P \V 3 €t( x > u ) 
where, as we have explained, 



oente[o,T] j=Q x&R d 



etc. H T is a Banach space under || • ||. 

In this section, if ( = (£*), where ( l G L 2 (Q, J 7 , P), then we define 
K(C) S = ^ P {CI^} and Z(C) = Dfl(y(C))- Therefore, for any t > 0, if 
C G L 2 (fi, Ji, P), then (T(C) S , Z(C)*) is the unique solution of the BSDE 

dY(() s = Z(() s .dB s , Y(C) t = (. 

It is easy to see that, if £ G //, . then Y(V k Q = V k Y(Q and Z(V k Q = 
V k Z(£ it ) for any t < T and k — 0, 1, 2. This follows from the fact that the 
mappings £ — > F(C) and £ — >■ Z(£) are both affine. 

Let £ G ff / . Then, according to the non-linear Cameron-Martin formula, 
for any fixed t < T, we define a probability Q 4 x on J^) by 



d Qt,x 



dP 



= exp 



/(y(e t (-,x)).).dB,-i / i./fn-^c.-nj./.s 
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We will omit the argument • (a sample point) and the space variable x if no 
confusion may arise. Under Q t x , B(£ t ) a = B s — J Q S f(Y(£ t ) r )dr (0 < s < t) 
is Brownian motion up to time t. Define 

X(0t = x + B(C t ) t 

= x + Bt- [ f(Y(Q s )ds for t < T. 



According to Theorem 13.11 we want to find a fixed point £ G H?'- £ 
Uq(X(£)). To this end we define $(£.) = Uq(X(£).) for any £ G // / Then 

w(0* = vu (x(0t)vx(0t 

and 

v 2 $(a = v 2 u (x(0t)(vx(0t,vx(0t) 
+vu (x(0t)v 2 x(0t- 

Lemma 5.1 l)For any T > and £ G i3Y 

WHO \\<C (l + d) 2 + OAT {(2d +!) + (! + CiT) | iei | } | iei i . 
2)IfK = 2C (l + d) 2 and 

T < 1 A 1, 



V d + 1/2 + C (l + Ci) (1 + df 
then ||$(0II <K as long as £ G // / and ||£|| < K. 

Proof. By definition for any t < T and any x (but the argument x is 
suppressed from the notations for simplicity, and | • |oo denotes the essential 
supremum norm) 

VX(0t = 4« - f t Vf(Y(z t ) a )Y(y£ t ) a ds 

Jo 

and 

-t 

72 



V 2 X(0 t = -/ ^f(Y(Q s )(Y(VQ s ,Y(VQ s )ds 
Jo 

- [ Vf(Y(Q s )Y(V%) s ds. 
Jo 
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From these equations and the fact that the conditional expectation is a con- 
traction on L°°(Q, T t -, P), one can easily to see the following estimates: 



t\ooi 



and 

Therefore 



and 



|V 2 X(O t | ^CxflV&f + CxflV 2 ^. 
1^(01 < Co, 

|v$(0*l < c |vx(o*| 

< C {d+C7i*|Ve t |oo} 



|V 2 $(0t| < Co|VX(O t | 2 + Co|V 2 X(O t | 



< Co {d + dt| V^loo} 2 + QAt [| + I V^loo] 

which yield the required estimates. ■ 

Lemma 5.2 Let T > 0. There is positive constant depending only on C , Ci 
and d sncn £na£ 

\\m - m\\ < Koani+T) (i + neii + ini) lie - «n 

/or any £,n G if T . 

Proof. By a simple computation, we have 



\x(0t-x( v ) t \ < 



f f(Y(Ct-Vt)s)ds 
Jo 



\VX(0t - VX(n) t | < 



+ 



(v/(y&).) - v/(y(n t ) s )) Y(vQ.d8 

*Vf(y(Tl t ) a )Y(V(t t -Tl t )).d8 

< d /V(e t -i7 t ) a ||y(ve t ).|ds 
+d | |y(v(e t -i? t )).|ds 
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and, similarly 

|V 2 X(0«-V 2 X(i7) t | < d^.-^iiv^r 

+Cl t(\V^\ + \V Vt \)\V^ t -V Vt \ 

+Cit\£ t - ^||v 2 e,| + c,t\ v% - v\| 

and the estimate follows easily from these inequalities. $(£.) = uq(X(£).) 

|^)t-^)t|<C C a ^ t -%U 



|V$(Ot-V$(r/) t | = C |X(Ot-X(^) t ||VX(^| 

+c |vx(e)t - VX{ V ) t \ 

< c^dti^-^u^+dtivetU) 
■K7 Cit(ie t -%iive t i + |v^-vi ?t i) 

■ 

Now we are in a position to establish a Bismut type formula (see Bismut 
PQ for the linear case) for the solution of quasi-linear system (I4.ip . 

Theorem 5.3 There is T > depending only on d, Cq and Ci, so that there 
is a unique fixed point £ of Q in Ht- Moreover 

V^(x,t) = E p (V^(-,x)) ; j = 0,1,2. (5.1) 

u(x + B(£ t )„, t-s) = Y(£ t ) s for all s <t <T a.e. 

and 

Vu{x + B(Q S , t-s) = Z(Q S for all s <t <T a.e. 

Proof. According to Theorem [3TT| for t < T we have u(x, t) = i? p (^(-, x)|J : o)• 
Taking expectation we obtain m(x, t) = E p (£ t (-,x)). Since £ t (u,-) G cf(R d ,R m ) 
so we may take derivatives under integration to obtain ( 15. ip . ■ 
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